PROCEEDINGS 

°F SCIENCE 



o 

D 
(N 



o 
u 

6 



in 
o 

O 

o 



X 



Preparation to the CMB PLANCK data analysis, 
estimation of tlie contamination due to the galactic 
polarized emissions 



L. Fauvet* 

LPSC, Universite Joseph Fourier Grenoble 1, CNRS/IN2P3, Institut National Polytechnique de 
Grenoble, 53 avenue des Martyrs, 
38026 Grenoble cedex, France 
E-mail: 



f auvetSlpsc .in2p3.fr 



J.-F. Maci'as Perez 

LPSC, Universite Joseph Fourier Grenoble 1, CNRS/IN2P3, Institut National Polytechnique de 
Grenoble, 53 avenue des Martyrs, 
38026 Grenoble cedex, France 
E-mail: }nacias@lpsc .in2p3.fr 



This work is point of the preparation to the analysis of the PLANCK satelhte data. The PLANCK 
satellite is an ESA mission which has been launched the 14th of may 2009 and is dedicaced 
to the measurement of the Cosmic Microwave Background (CMB) in temperature and polariza- 
tion. The presence of diffuse Galactic polarized emissions disturb the measurement of the CMB 
anisotropics, in particular in polarization. Therefore a precise knowledge of these emissions is 
needed to obtain the level of accuracy required for PLANCK. In this context, we have developed 
and implemented a coherent 3D model of the two mains polarized Galactic emissions : syn- 
chrotron and thermal dust. We have compared these models to preexisting data: the 23 GHz band 
of the WMAP data, the 353 GHz Archeops data and the 408 MHz all-sky continuum survey. We 
extrapolate these models to the frequencies where the CMB dominates and we are able to estimate 
the contribution of polarized foreground emissions to the polarized CMB emission measured with 
PLANCK. 
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1. Introduction 

The PLANCK satellite, currently in flight, should give the more accurate measurement of 
the anisotropics of the CMB in temperature and polarization with a sensitivity of IjjiK and an 
angular resolution of 5 arcmin [||]. In particular its estimation of the BB -modes should set an 
upper limit on the tensor-scalar ratio (expected at 0.1 [^). The knowledge of this ratio should 
confirm the existence of primordial gravitational waves generated during the inflation and would 
set the energy scale of the inflation [ p3| ] and provide constraint on inflationnary models [[T]]. In 
order to obtain its optimal sensitivity it is required to estimate the foreground emissions and the 
residual contamination due to these foreground emissions on the CMB signal. Indeed for the full 
sky these emissions have the same order of magnitude than the CMB in temperature and dominate 
by a factor of 10 in polarization [^. The principal polarized Galactic microwave emissions come 
from 2 effects : thermal dust emission and synchrotron emission. The synchrotron has already 



been measured by the 408 MHz all-sky continuum survey 1 1 1], by Leider between 408 MHz and 



1.4 GHz [25], by Parkes at 2.4 GHz M, by the MGLS Medium Galactic Latitude Survey at 1.4 



GHz | |24| ] and by the satellite WMAP Wilkinson Microwave Anisotropics Probe (see e.g. []12[]). 
The synchrotron emission is due to ultrarelativist electrons spiraling in a large-scale magnetic field 
and is dominant at low frequencies. The dust thermal emission which have already been well 
constrained by IRAS [^], COBE-FIRAS ^ and Archeops [|T6l, g] is due to dust grains which 
interact with the Galactic magnetic field and emit a polarized submillimetric radiation [Q] and 
dominates at high frequencies. The polarization of these two radiation is orthogonal to the field 
lines. To obtain a realistic model of these emissions we propose models based on a 3D modelling 
of the Galactic magnetic field and of the matter density in the Galaxy. The models are optimized 
using preexisting data and then are used to estimate the bias due to these emissions on the CMB 
measurement. 



2. 3d modelling of the Galaxy 



A polarized emission is described by the Stokes parameters I, Q and U [jlij]. For the polarized 



foreground emissions integrating along the line of sight we obtain, for synchrotron D20| ] : 



u. 
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(2.1) 
(2.2) 
(2.3) 



where B„, Bi and B, are the magnetic field components along, longitudinal and transverse to the 



ligne of sight, is the polarization fraction set to 75% [20]. Inas is a template temperature map 



obtained from the 408 MHz all-sky continuum survey [11]. The maps are extrapolated to the Planck 
frequencies using the spectral index jS^ which is a free parameter of the model. 
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For the thermal dust emission : 



[353) J"" n, ' ('-^^ 

\353) J"' ' ^^-^^ 

where the polarization fraction pd is set to 10 % [|l9|], j3j is the spectral index (set at 2.0) and fnorm 
is an empiric factor, fit to the Archeops data. The Igfd map is the model 8 of [|8|]. 



Id 


= I^fd 


Qd 


= Isfd 


Ud 


= hfd 



The models are based on a exponential distribution of relativistic electrons on the Galactic disk 
following [||] where the radial scale /z,- is a free parameter. The distribution of dust grains na is also 
choose exponential [^. The Galactic magnetic field is composed of two parts: a regular component 
and a turbulent component. The regular component is based on the WMAP team model [18] which 
is close to a logarithmic spiral to reproduce the shape of the spiral arms The pitch angle p 

between two arms is a free parameter of the model. The turbulent component is described by a law 
of Kolmogorov [Q, spectrum of relative amplitude Afurh- 



3. Comparison to data 



We computed Galactic profiles in temperature and polarization for various bands of longitude 
and latitude and various values of the free parameters. In order to optimize these 3D models 
we compare them to Galactic profiles computed from preexisting data using a test. For the 



synchrotron emission in temperature, we use the 408 MHz all-sky continuum survey [ ]11[ ] as shown 
on Figure [T]. In polarization we compared to the K-band WMAP 5 years data. Thermal dust 
emission model is optimized using the polarized Archeops data [§] at 353 GHz. 

The best fit parameters for the 3D model in polarization are given in the Table |l|. The results 
are consistent for the 3 sets of data, and in particular we obtain compatible results for the syn- 
chrotron and thermal dust emission models. Aturb is not strongly constrained but its range of best 
fit value is compatible with previous results |23, |l^, 10 1. hr is badly constrained as was already the 



case in Sun etal[13\. The best fit value of the pitch angle p is compatible with results obtained by 



other study [ p3| , [Tq ]. The best fit value for the spectral index of the synchrotron emission is lower 
than value found by [| 



ISp but it is probably due to the choice of normalisation using the 408 



MHz template. 



Using the best fit parameters obtained for the Galactic emissions models we computed maps 
and power spectra in temperature and polarization for synchrotron and dust thermal emission. We 
compare them to maps and power spectra built respectively using polarized WMAP and ARCHEOPS 
data. Like represented in the Figure ^ the synchrotron emission model is efficient to reproduce the 
global feature of the data in polarization. We show on Figure ^ the angular power spectra computed 
from the Archeops data at 353 GHz and the thermal dust model using the method presented in [p^]. 
Our model efficient to reproduce the features of the spectra at all scales. 
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Figure 1: Galactic profiles in temperature at 408 MHz built using Haslam data in black and our synchrotron 
emission model for various values of the pitch angle p (form green to red). 
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Figure 2: Clockwise from top left : power spectra Cf^, Cf^, Cf^, Cj^ , Cj^, Cf^ at 23 GHz built with the 
WMAP 5 years data (black) and the model of synchrotron emission with BSS magnetic field for the best fit 
model parameters (green to red for a spectral index between -3.2 and -3.4), applying a Galactic cut \b\ < 5°. 
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Table 1: Best fit parameters for synchrotron and thermal dust emission models and 3 a confidence levels for 
the best fitting model. 
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Figure 3: Clockwise from top left : power spectra C/ ' C, , C,, C,, at 353 GHz computed 

from Archeops data (black) and the model of thermal dust emission with BSS magnetic field for the best fit 
model parameters f red) without applying Galactic cut. 

From the above best fit parameters we estimate the contamination of the CMB PLANCK data 
by the polarized galactic emissions. Figure ^ shows the temperature and polarization power spec- 
tra at 143 GHz for the CMB ' (red) and the Galactic foreground emissions applying a Galactic cut 
of 1^7 1 < 15°. The residual foreground contamination seems to be weak but for the BB-modes for 
which an accurate foreground substraction is extremely important for the detection of the primor- 
dial gravitational waves. 

4. Conclusions 

We propose in this study consistent models of the main Galactic polarized emissions based 
on a 3D modelisation of the Galaxy. By comparison with preexisting data we are able to give 

' We simulate CMB assuming cosmological parameters for a model ACDM like proposed in | [l3| ] with a ratio tensor- 
scalar of 0.03. 
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Figure 4: Clockwise from top left : power spectra Cf^, Cf^, Cf^, Cf^, C™, Cf ^ at 143 GHz for \b\ < 15° 
(see text for details). 



consistent constraints on the parameters of the synchrotron and dust thermal emissions models 
compatibles with thus appear in the literature. From this we build map and power spectra enable to 
reproduce the features of the data at various frequencies. Using a rough mask we then estimate the 
residual contamination due to these foregrounds on the assumed CMB PLANCK data. 
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